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We show that at one-loop order, negative-dimensional, Mellin-Barnes' (MB) and Feynman 
parametrization (FP) approaches to Feynman loop integrals calculations are equivalent. Start- 
ing with a generating functional, for two and then for iV-point scalar integrals we show how to 
reobtain MB results, using negative-dimensional and FP techniques. The N— point result is valid 
for different masses, arbitrary exponents of propagators and dimension. 
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^ ' I. INTRODUCTION 
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The amazing comparison[l] between experimental determination and theoretical prediction of the anomalous mag- 
netic momentum of the electron, is the greatest motivation - in our opinion - to study and develop techniques that 
allow the precise calculation of higher order Feynman loop integrals. Recently there are also interest in studying 
00 ■ process like [2] e + e~ — > 3 jets and e + e~ — > 4 jets , so loop integrals with five and six external legs must be known. 

Physicists' battle against the tricky Feynman loop integrals is fought in many fronts. We can cite some of them: 
integration by parts method[3] seems to be the most powerful one, since one can in most cases reduce the number 
of loops, e.g., the scalar massless two- loop master can be rewritten as a sum of two simpler integrals: a two- loop 
self-energy with an insertion[4, 5] plus the square of a one-loop self-energy. This is a very simple example - where one 
must deal with a greater number of simpler integrals with powers of propagators shifted - of a powerful technique, 
see for instance a 5-loop calculation in [6]. 

Integration by parts is used also associated with another methods. In fact, one can not evaluate a Feynman loop 
f"*) |' integral using the above mentioned technique alone. It simplifies the original diagram but does not solve it. In 
order to carry the integration out Gehrmann and Remiddi [7] did use the differential equation method, introduced 
by Kotikov[8], and solved a large class of difficult of problems. Glover and collaborators completed the study of the 
whole class 2 — ► 2 of two-loop scattering [9]. Also, Gegenbauer polynomial method has been used in order to study 
complicated process[10]. 

Other methods that make use of decomposition of complicated integrals, like the one-loop pentagon[ll], were 
developed as well as string inspired ones [12]. See also [13-16] for other important approaches, a very powerful 
numerical technique on [17] and for a review on the progress of loop calculations [18]. 

Mcllin-Barnes(MB) and negative dimensional integration method (NDIM) are two other interesting and powerful 
techniques to tackle such Feynman integrals. Mellin-Barnes approach relies on the relation, 
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in other words, we rewrite each propagator as a Mellin transform. However, these parametric integrals are not 
difficult to solve - as it happen in the Feynman parametrization approach where the integrals become more and more 
complicated - because one can apply Cauchy theorem and two Barnes' lemmas [19]. The MB approach is being greatly 
used by Tausk[20], Smirnov[21], Davydychev[22, 23] and co-workers in order to tackle two and three-loop integrals. 
The results are always expressed as generalized hypergeometric functions which depend on adimensional ratios of 
momenta and/or masses, space-time dimension D and exponents of propagators. 
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On the other hand NDIM is a technique whereby it is not necessary to introduce parametric integrals, the Feynman 
integral is the result between the comparison of two calculations: a gaussian-like integral (the generating functional of 
negative-dimensional integrals) and a Taylor expansion of the generating functional[24]. It is worth mentioning that 
in the NDIM context performing the calculation for a particular set of exponents of propagators present the same 
difficulties than perform the same task for arbitrary values of them. The results are, just like in the MB approach, 
given in terms of generalized hypcrgcomctric functions which depend on the same quantities mentioned above. 

One can then rightfully ask: is there any connection between these two approaches? The answer is yes. The propose 
of our paper is to show that they are equivalent, at least at one-loop order. In fact, one could argument that the 
results must be the same if one correctly applies both methods, we will show this explicitly. However, it can be useful 
and interesting to study which of them is more powerful when the number of external legs increase. In one hand 
NDIM demands computer facilities in order to solve a large number of systems and browse the big number of results; 
on the other hand MB does not require computers but the integrals must be calculated one by one. Another point to 
observe is that NDIM relics on grassmannian integrals and MB on Mellin transform, i.e., apparently disconnect (as 
far as we know) subjects. Also, showing the equivalence between them and knowing the routes which take to NDIM 
or to MB one could build, for instance, a technique like NDIM in order to tackle problems in finite temperature field 
theory, like calculations of heat-kernel which can be dealt with using Mellin integrals. 

The outline for our paper is as follows: in section II we present a step-by-step calculation of 2-point scalar integral 
starting with NDIM approach and arriving at an expression originally obtained by Davydychev[23] using MB approach, 
then we repeat the same process using FP. In the next section, we deal with an arbitrary number N of external legs, also 
starting in the negative-dimensional approach and showing how to obtain the Davydychev's original result calculated 
in the MB scheme; we carry out the same the integrals with FP. In the final section we present our conclusions and 
a discussion concerning the three methods. 



II. ONE-LOOP 2-POINT FUNCTION 



In this section we present the calculations to evaluate the one-loop two-point scalar integral within NDIM scheme 
and compare this result with the obtained by Mellin-Barncs approach. Consider the integral, 
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which is the usual generating functional for negative-dimensional integrals. We will always begin with this kind of 
generating functionals, for two and for n-point scalar integrals, and after some manipulations arrive at results which 
were obtained previously by other authors, using MB approach. 

The first step in NDIM context is a series expansion of the above integral, 
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and then a multinomial expansion gives us, 
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where (a)& = T(a + b)/F(a), is the Pochhammer symbol. 

Compare the power of the parameters a and (3 between (3) and (9) we have the following constraint equations, 



that after solving for jo and ji we have 



ai = -D/2 + jo -bi- ji 
a 2 = h + ji, 



jo = ai + 0.2 + D/2 = (T 2 
ji = ai - h. 



Performing the substitution of this result in (9) and the analytic continuation to a\,a 2 < 0, we arrive at, 
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which is the well-known result for 2-point scalar integrals, with different masses, see [19]. 



A. Two-point function via Feynman Parametrization 



The most popular technique to deal with loop integrals is certainly, Feynman parametrization. It is the one the 
students learn on field theory courses, and one of the few textbooks introduce (the other is a— parametrization). 
Depending on the manipulations one performs, it can turn the original loop integrals into a hefty one. We will 
proceed in a slightly different route. Our aim is to show how one can obtain the previous results for 2-point functions, 
given in terms of hypergeometric functions, using FP since in most cases the results calculated through FP are written 
as poly logarithms, Li n (z), n = 0, 1, 2, 3, 4. 

Hypergeometric functions have an advantage over dilogarithms, for instance, in the case of photon-photon scattering 
scalar integrals. The result for |s/4to 2 < 1, \t/4m 2 \ < 1 can be written as a single Appel function F3 of two variables 
and 5 parameters, on the other hand, the same result can be recast as a sum of several Li 2 {zj) functions of complicated 
arguments Zj, see for instance [25]. 

Consider the function 
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where 0,1,0,2 > 0. We note that when xq = 1, x 2 = we have the well-known Feynman paramctrization to the 
propagator of , that is = . Such modification will turn simpler the calculation of N— point integrals in 
section IIIA. This expression can be rewritten as follows, 
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that after the evaluation of the integral in k using the well-known formula, 
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These integrals in x\ can be evaluated, with xq = 1, x 2 = 0, using 
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which is exactly the former result (14). This result show that, to one-loop two-point, the NDIM, Feynman parametriza- 
tion and Mcllin-Barncs representation are equivalent. The other kinematical regions can be obtained through analytic 
continuation of hypergeometric function above (see [19, 22, 26]). 



III. N-POINT FUNCTION 



In this section we present the generalization of the previous ideas in order to the obtain the Mellin-Barnes result 
for the scalar integral associated to n-point function. We consider a one-loop Feynman diagram with n external legs 
with momenta: p\ = l 2 — h,P2 = h — h, ■■■,Pn = h — l n , and internal momenta k — l\,k — l 2 , ...,k — l n . From a similar 
reasoning, we begin with the generating functional, 
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where (li, mi, a\\ Z n , , m n , a n ) represents the n— point functions for negative values of ai and is given by 
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where ijj ~ h — lj ■ After a new expansion in the right side of the expression above, we have 
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We need to do now the comparison term-by-term between «j powers in eq.(22) with the ones of eq.(19). We obtain 
fi = ^ and the solution of system above will be given by, 
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Performing the substitution of the solutions above in (22), we arrive at 
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that after carrying out analytic continuation to negative values of the a, provides, 
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The above result in also an expression for the n-point scalar integrals with arbitrary exponents of propagators and 
dimension, in the MB scheme. However, it was not this one the obtained by Davydychev in [23]. Formula (27) is a 
new result. 

It is important to observe that the above result, eq.(27), is valid since the series is convergent which means \l1j/l\ 2 \ < 
1 and |m?/Zi 2 | < 1 5 that is, external momentum greater than masses. Conversely, Davydychev's result holds in 
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another kinematical region, namely, where llfj/m^l < 1 and |1 — mf/m n \ < 1, i.e., when masses are greater than 
incoming/outcoming momenta. This result could be obtained, in principle, from the Davydychev's formula through 
analytic continuation. However we stress the point that such analytic continuation formulas are not known for multiple 
hypergeometric series (in general these formulas are known only in the case of single and double series). 
Other form to represent the n-point function can be obtained also from expansion of (21), that is 
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Analytically continuing and performing the substitution of the result above in (28) and compare its on powers with 
(19), we arrive at, 
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that after the analytic continuation to negative values of a i7 we get 
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which was the result obtained by Davydychev in [23] using MB approach. So, with the generating functional (18) as 
the starting point — the same which we have, in a previous paper [27] in the NDIM approach, used to show how to 
obtain a general formula to any scalar one-loop Feynman integrals, in covariant gauges — we were able to reproduce 
a MB result eq.(33) and more, to present another formula (27), also valid for n-point scalar integrals. 
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A. N-Point function via Feynman Parametrizat 

Our final task in this paper is to show how to solve an TV-point scalar integral using FP technique. As far 
know there is no such result in the literature calculated using FP. Of course, it has to be the same we obtained 
using NDIM and Davydychev's [23] MB approaches. 
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where a, > 0. This function for xo = 1 i n = 0, represent the Feynman parametrization to the integral of type (20). 
The integral in k above can be evaluated using (15), 
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Using now 
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performing also the Taylor expansion of the argument of the integral above, we get 
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The integral above can be evaluated with help of (16). If we take xq = 1, x n = 0, — J^ n \ we arrive 
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This result is the same one obtained in the previous subsection via NDIM in (33). This agreement show that, in 
one-loop level, the NDIM, Feynman paramctrization and Mcllin-Barnes representation present the same results and 
are equivalent: all of them can be used to solve all scalar Feynman loop integrals at one-loop order, with general 
masses, arbitrary exponents of propagators and dimension. 



IV. DISCUSSION AND CONCLUSION 



So far we have made calculations in order to show that the same class of generating functionals can be used to 
reproduce MB results. Depending on which Taylor expansions one carries out one can proceed in the NDIM or MB 
routes. The final results will be, obviously, the same, given in terms of generalized hypergeometric functions, being 
the exponents of propagators and space-time dimension arbitraries. 

However, one could ask which of these two routes, if any, is the one where Feynman integrals become simpler to solve. 
The first point to observe what are the tools one has to master in order to tackle such integrals in both approaches: 
contour integration, Cauchy theorem and Barnes' lemmas for MB, and solving system of algebraic equations for 
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NDIM. So far, so good. Second, the results, despite they will be the same, have to be worked out one-by-one in 
the MB context, on the other hand, using NDIM and solving the system of algebraic equations gives one all the 
possible solutions (generalized hypergeomctric functions) for the Feynman integral in question. Group them is a 
straightforward task: linear independent functions have to be summed, each set is a possible result in a given region 
if convergence [26]. Third, the massless case needs to be known in the case of MB in order to tackle massive integrals; 
not so in NDIM. 

We can summarize both approaches in following table, 
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Result: system of algebraic equations 
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Use the results of the systems 


10 


and Barnes' lemmas 


Analytically continue to positive D 


11 


One have one final result among several 


One have all the scries (final results) 



in the step number 3 one can proceed as we have done in this paper, expanding the exponential, or as we did in 
our previous works taking a Taylor expansion for each argument of the exponential. The final step, 11, is to be 
understood in the following manner: in order to obtain all possible generalized hypergeomctric functions (which come 
in NDIM) using MB one has to repeat the above procedure choosing other sequence of contours, we mean for instance 
left-left-right-left-right and another one left-left-right-right-right, these two can give, in principle distinct generalized 
hypergeometric series. Some of them, will of course result in zero, since there can be no poles inside the contours. 
These ones are also contained in the NDIM approach, since some determinants can vanish, a much simpler calculation 
that can be implemented in softwares like Mathematica. 

The textbook technique, FP, can be made simpler if one introduces two extra parameters xo and x n , and takes 
series expansions in the parameters (xq — xi), (x\ — X2), (x n -i — x n ). In the end of the day one makes xq = 0, 
x n = 1 and uses the well-known beta function integral representation. Then, the remaining expression is the result 
written as a generalized hypergeometric function. 

A. Conclusion 

We have shown that negative-dimensional integration method (NDIM), Feynman parametrization (FP) and Mcllin- 
Barnes' approach to scalar Feynman loop integrals, at one-loop level, give the same results. It depends only on how 
one choose to Taylor expand the generating functional (18). We present detailed calculations for two-point scalar 
integrals, with arbitrary masses, exponents of propagators and space-time dimension (in covariant gauges). Then we 
tackle a general scalar iV-point integral, with different masses, and did show that the general formula of Davydychev 
[23] and ours [27] agree, as well as another one obtained via FP worked out, as far as we know, for the first time. It 
is our opinion however, that NDIM is simpler than MB, since all the possible results for the integral in question are 
obtained simultaneously, and in MB they must be calculated one by one, or through analytic continuation formulas, 
if such formulas were known, depending on the hypergeometric functions. FP is also a very powerful technique if one 
introduces two extra parameters and take Taylor expansions properly. In doing so, FP can become even simpler than 
NDIM, since one obtain the full result and does not have the drawback of searching among a huge amount of possible 
solutions. 
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